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Abstract 

This paper develops a rigorous asymptotic expansion method with its numerical 
scheme for the Cauchy-Dirichlet problem in second order parabolic partial differential 
equations (PDEs). As an application, we propose a new approximation formula for pric- 
ing a barrier option under a certain type of stochastic volatility model including the 
log-normal SABR model. Q 
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1 Introduction 

Numerical methods for the Cauchy-Dirichlet problem have been a topic of great interest 
in stochastic analysis and its applications. For example, in mathematical finance the Cauchy- 
Dirichlet problem naturally arises in valuation of continuously monitoring barrier options: 

Cb™(T,x) = E[f(X§.)l {r>T} ] = E[f(X*)l {miBte[QT]Xt>L} }. (1.1) 

Here, T > is a maturity of the option, and (Xf) t denotes a price process of the underlying 
asset starting from x (usually given as the solution of a certain stochastic differential equation 
(SDE)). Also, L stands for a constant lower barrier, that is L < x, and r is the hitting time 
to L: 

T = inf{t G [0,T] : X x t < L}. (1.2) 
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It is well-known that a possible approach in computation of CBa,rrier(T,x) is the Euler scheme, 
which stores the sample paths of the process (X*)t through an n-time discretization with the 
step size T/n. In applying this scheme to pricing a continuously monitoring barrier option, 
one kills the simulated process (say, (Xf.)i) if X*. exits from the domain (L, oo) until the 
maturity T. The usual Eular scheme is suboptimal since it does not control the diffusion 
paths between two successive dates tj and t i+ i. the diffusion paths could have crossed the 
barriers and come back to the domain without being detected. It is also known that the error 
between CBarrier(^, x) and C , Barrier(^ 1 , x) (the barrier option price obtained by the Euler scheme) 
is of order yT/n, as opposed to the order T/n for standard plain-vanilla options. (See [7J) 
Therefore, to improve the order of the error, many schemes for the Monte-Carlo method have 
been proposed. (See [16] for instance.) 

One of the other tractable approaches for calculating CBamer{T,x) is to derive an analyt- 
ical approximation. If we obtain a closed form approximation formula, then it is a powerful 
tool for evaluation of continuously monitoring barrier options because we do not have to rely 
on Monte-Carlo simulations anymore. However, from a mathematical viewpoint, deriving an 
approximation formula by applying stochastic analysis is not an easy task since the Malli- 
avin calculus cannot be directly applied, which is due to the non-existence of the Malliavin 
derivative D t r (see [I]) and to the fact that the minimum (maximum) process of the Brownian 
motion has only first-order differentiability in the Malliavin sense. Thus, neither approach in 
[TT] nor in [T^] can be applied directly to valuation of continuously monitoring barrier options 
while they are applicable to pricing discrete barrier options. (See [18] for the detail.) 

In this paper, we propose a new general method for approximating the solution to the 
Cauchy-Dirichlet problem. Roughly speaking, our objective is to pricing barrier options when 
the dynamics of the underlying asset price is described by the following perturbed SDE: 

/ dXr = b(X^, e)dt + a(X £ >*, e)dB u . . 

\Xr = x, 

where e is a small parameter, which will be defined precisely later in the paper. In this case, 
the barrier option price (11 .11) is characterized as a solution of the Cauchy-Dirichlet problem: 

d 

— u £ (t, x) + Jzf e w e (t, x) = 0, (t, x) G [0, T) x (L, oo), 

u £ {T,x) = f(x), x>L, ( L4 ) 

u £ (t,L) = 0, te[0,T], 

where the differential operator «5f e is determined by the diffusion coefficients b and a. Next, 
we introduce an asymptotic expansion formula: 

u £ (t, x) = u°(t, x) + ev^t, x) + ■ ■ ■ + e n-1 vJi-i(t, x) + 0(e n ), (1.5) 

where O denotes the Landau symbol. The function u°(t, x) is the solution of fll.4p with e = 0: 
if b(t, x, 0) and a(t, x, 0) have some simple forms such as constants (as in the Black-Scholes 
model), we already know the closed form of u°(t, x) and hence obtain the price. Then, we are 
able to get the approximate value for u £ (t, x) through evaluation of v®(t, x), . . . , x). In 

fact, they are also characterized as the solution of a certain PDE with the Dirichlet condition. 
By formal asymptotic expansions, (11.51) as well as 

J^ £ = + + ■ ■ • + e"" 1 ^! + ■ ■ • , 
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we can derive the PDEs corresponding to v k (t,x) of the form: 

JU°(t, x) + Sf°v° k (t, x) + g° k (t, x) = 0, (t, x) G [0, T) x (L, oo), 

u°(7» = 0, x>L, (L6) 

vg(*,L)=0, te[o,T], 

where g k (t, x) will be given explicitly later in this paper. Moreover, by applying the Feynman- 
Kac approach, we are able to obtain their stochastic representations. We will justify the above 
argument in a mathematically rigorous way with necessary assumptions in Section [5J 

The theory of the Cauchy-Dirichlet problem for this kind of second order parabolic PDE 
is well understood in the case of bounded domains (see [5], [B] and [2] for instance). As for 
an unbounded domain case such as (ll.4p . [TT] provides the existence and uniqueness results 
for a solution of the PDE and the Feynman-Kac type formula (cited as Theorem [1] below) . 
However, some mathematical difficulty exists for applying the results of [17] to the PDE ( II. 5J) . 
More precisely, the function g k (t,x) may be divergent at t — T. (If g k (t,x) is continuous on 
[0, T] x [L, oo), the existence and uniqueness of (jl.fip are guaranteed: see [5].) To overcome this 
difficulty, we generalize the Levi's parametrix method (which is used to construct a classical 
solution of the PDE) in Theorem [2j Furthermore, we show another representation of v k (t,x) 
by using the corresponding semi-group in Section [31 We notice that such a form is convenient 
for evaluation of v k (t,x) in concrete examples. 

In Section m we apply our method to pricing a barrier option in a stochastic volatility 
model: 

dS £ t = (c - q)S £ t dt + a £ t S £ t dBl S £ = S> 0, 



da £ t = e\(9 - a £ t )dt + evo*(pdB\ + y/l - p 2 dB 2 t ), a £ = a > 0, 

where c, q > 0, e G [0,1), X,9,u > 0, p G [—1,1] and B = (B X ,B 2 ) is a two dimensional 
Brownian motion. Then, we obtain a new approximation formula: for x = log S, 

^Barrier ( T ' e2 ) = E [/(St*) l{min < t <T S$>L}] 

~ P°f(x)+e [ P^_ r ^ P r D f(x)dr, 
Jo 

where Sf' x = exp(X t e,:c ), (Pf)t is a semi-group defined in Section [31 / is a payoff function 
and f(x) = f(e x ). Here, P£f{x) is regarded as the down-and-out barrier option price, 
C£^ rier (T, e x ) in the Black-Scholes model. Moreover, we confirm practical validity of our 
method through a numerical example given in Section [U Notice also that our example does 
not satisfy the assumptions introduced in Section [21 Thus, we generalize our main result 
and present weaker (but a little bit complicated) version of the assumptions in Section 15.11 
Furthermore, Section I5.2H5.4I list the proofs of our results. 

Finally, we remark that in the contrast to the previous works ([2], [3], [8], [9] for ex- 
ample), which start with some specific models (the Black-Scholes model or some type of fast 
mean-reversion model) and derive approximation formulas for (discretely or continuously mon- 
itoring) barrier option prices, we firstly develop a general asymptotic expansion scheme for the 
Cauchy-Dirichlet problem under multi-dimensional diffusion setting; then, as an application, 
we provide a new approximation formula under a certain class of stochastic volatility model 
that can be widely applied in practice (e.g. in currency option markets). 
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2 Main Results 



Let b : R d x I — ► R d and a : R d x / — ► R d ®R m be Borel measurable functions (d, m E N,) 
where / is an interval on M including the origin (for instance / = (—1, 1).) We consider 
the SDE ( II. 3p for any i6R d and e E I; we will introduce the assumptions for existence and 
uniqueness of a weak solution of (jl.3p later. 

We are interested in evaluation of the following: for a small e, 



u e (t, x) 

T-t 

(t, x) E [0, T] x D (2.1) 



E 



exp (- c(X £ r ' x , e)dr^j f(X^ t )l{r D (x^)>T-t} 



for Borel measurable functions / : R d — > R and c : R d x I — > R, a positive real number 
T > and a domain D C R d ; D C R d is the closure of D and r D (w), w E C([0,T]; R d ), stands 
for the first exit time from D, that is 

t d (w) = inf{t E [0,T};w(t) (£ D}. 

As mentioned in Section the right-hand side of (12.11) corresponds to a barrier option price of 
knock-out type with maturity T in finance. We regard (X^' x ) t as the underlying aseet prices 
and the expectation E[-] is taken under a risk-neutral probability measure. The boundary dD 
of the domain means the trigger points of the option and / represents a payoff at maturity. 
The function c represents a short-term interest rate. Our setting includes the case of D = R d , 
which corresponds to a price of an European option: 



u £ (t,x) = E 



exp t c(Xp x ,e)dr\ f{X e T 



£,X > 
tJ 



For applications to option pricing, see Section H] for the details. 
Now we introduce our assumptions. 

[A] There is a positive constant A\ such that 

\a ij {x } e)\ 2 + \b\x,e)\ 2 < A 1 (l + \x\ 2 ), x E R d , eel, i,j = l,...,d. 

Moreover, for each e e I it holds that cr JJ, (-, e), b l (-,e) e L for i, j — 1, . . . , d, where C is 
the set of locally Lipschitz continuous functions defined on R d . 

[B] The function f(x) is continuous on D and there are Cf > and m E N such that 
|/(aOI < C f0- + k! 2m ) ; x e K d . Moreover, f(x) = on R d \ D. 

Note that under [A], the existence and uniqueness of a solution of (11.31) are guaranteed on 
any filtered probability space equipped with a standard c?-dimensional Brownian motion, and 
Corollary 2.5.12 in [TU] and Lemma 3.2.6 in [13] imply 

E[ sup \X £ ' X - x\ 21 } < Cit l ~\l + |x| 2 '), (t,x) E [0,T] x R d , I e N (2.2) 

0<r<t 

for some C; > which depends only on / and Ai. Moreover, {X x ) r has the strong Markov 
property. Although the assumptions [A]-[B] are not always satisfied in our example in Section 
IH we can weaken them, and will introduce more general conditions in Section 15.11 
We continue to state our assumptions. 
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[C] There is a positive constant A 2 such that c(x,e) > —A 2 for x G D, e G I. Moreover, for 
each e G I, it holds that c(-,e) G C 

[D] The boundary <9-D has the outside strong sphere property, that is, for each x G dD there 
is a closed ball E such that E (1 D = <p and E (1 D = {x}. 

[E] The matrix (a u (x, e))ij is elliptic in the sense that for each e E I and compact set K C M. d 

d 

there is a positive number \i e K such that a^(x, e)C^ > /^e^l^l 2 f° r an Y x £ K and 

i,i=i 

£ G R d . In the case of e = 0, we further assume 

/ioiei 2 < J2 aij ( x Mv < /io'iei 2 , x g e g M d 

for some /i > 0. 
Let us define a second order differential operator Jzf £ by 



1 ^ 

^ £ = - V a* (x, e) 7: — — — : + (x, e) 



d 2 A,.-, , 

: — ClX, £ 



2 ^— ' dx l dx^ dx 

i,j=l i=l 



where a u = a lk a^ k . We consider the following Cauchy-Dirichlet problem for a PDE of 

k=l 

parabolic type 

d 

— u £ {t, x) + ££ £ u £ (t, x) = 0, (t, x) G [0, T) x D, 

u £ {T,x) = f(x), xeD, ( 2 - 3 ) 

u £ (t,x)=0, (t, x) G [0, T] x <9-D. 

The following is obtained by Theorem 3.1 in |17j . 

Theorem 1. Assume [A]-[E]. For each e G I, u £ (t,x) is a (classical) solution of (I2.3P and 

sup \u £ (t,x)\/(l + |x| 2m ) < oo. (2.4) 

(t,y)e[0,T]xD 

Moreover, ifw £ (t,x) is also a solution of (12. 3p satisfying the growth condition 

sup K(t,a;)|/(1 + |x| 2m ') < oo 

(t,2/)e[0,T]xD 

/or some m' G N, then u £ = w £ . 

To study an asymptotic expansion of u £ (t,x), we assume 

[F] For each i,j = l,...,d the functions a v (x,0), b l (x,0) and c(t,x, 0) are bounded on 
[0,T] x D, and there exist constants A3 > and a G (0, 1] such that 

\a^(x, 0) - <7*% 0)| + |6*(x, 0) - 6%, 0)| + |c(x, 0) - c(y, 0)| < A 3 \x - y\ a , x,y <E D. 



[G] Let n G N. The functions a JJ (x, s), b\x,e) and c(x,e) are n-times continuously dif- 
ferentiable in e. Furthermore, each of derivatives d k a l; ' /de k , d b l /de , d k c/de h , k = 
1, . . . , n — 1, has a polynomial growth rate in x 6 M 6 * uniformly in e G /. 

By [G], we can define Jzf fc °, /c G N, as 

~ n 1 fl^#d ,J , „, d 2 J^8 k b\ n . d d k c, n A 

fe k\ ) 2 ^ de k dx l dx? ^ de k y ' dx l de k K M v 7 

Our purpose is to present an asymptotic expansion such that 

u £ (t,x) =u (t,x)+ev° 1 (t,x) + --- + e n - 1 v°_ 1 (t,x) + O(e n ), e -»■ 0. (2.6) 
Here, v°(t, x), A; = 1, . . . , n — 1, are given as the solution of 
d 



^v&t, x) + Sf°v° k {t, x) + a°(*, x) = 0, (t, x) G [0, T) x D, 

^(T,x) = 0, xeD, 

v° k (t,x) = 0, (t,x) G [0,T] x dD, 



(2.7) 



where g"(t, x) is given inductively by 



fc-i 



S°(t, x) = J2f«°(t, x) + *2-iV?(t, *) ■ (2-8) 



i=i 

To state the existence of such a function v k (t, x), we prepare the set % m,a ' p of g G C([0, T) x 
I)) satisfying the following conditions: 

• There is some M 9 G C([0,T)) n L p ([0,T),dt) such that 

|g(t,x)|<M 9 (t)(l + |x| 2m ), £ G [0, T), x, ?/ G D. (2.9) 

• For any compact set K C D there is some M 9 ' K G C([0, T)) n L p ([0, T), dt) such that 

\g(t,x) - g(t,y)\< M 9 ' K (t)\x -y\ a , te[0,T), x,y G K. 

Then, we have the next theorem of which proof is given in Section 15.21 
Theorem 2. Assume [A]-[G\. Let g G 1-i m > a ' p for some p > 1/a. Then, the following PDE 
d 

—v(t, x) + JSf°u(t, x) + g(t, x) = 0, (t, x) G [0, T) x D, 

u(T,x)=0, x G -D, ^ 2 ' 10 ^ 

w(t,x) = 0, (t, x) G [0, T] x dD 

has a classical solution v such that 

\v(t,x)\ < C(l+ |x| 2m ) (2.11) 

/or some C > which depends only on a(-, 0), &(•, 0), c(-, 0), -D and M 9 . Moreover, if w is 
another classical solution of (I2.10p which satisfies \w(t,x)\ < C exp(/3|x| 2 ), (t,x) G [0, T] x5, 
/or some C', /3 > 0, £/ien v = w. 
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We also put the next assumption: 
[H] u° E g m > a ' p for some p > 1/a, where 



■>m,a,p 



{geC l > 2 ([0,T)xD)nC([0,T]xD) 



dg 

dx l ' dx l dxi 



hJ 



It is easy to see that the assumptions [F]-[H] imply g® E % mi > a ' p for some mi E N. 
Therefore (12. 7p with k = 1 has a unique classical solution under [A]-[H]. Similarly, if 
v®, . . . ,vl exist and are subject to g mk > a > p for some m k E N, then the unique classical solution 
f £ +1 of (12. 7p exists. We introduce our final assumption. 

[I] It holds that v° k E Q m ^ a 'P, k = 1, . . . , n - 1 for some m n e N. 

We remark that x) has the stochastic representation: 

*{T-t)Ar D (X ' x ) 



v° k (t,x) = E 



exp - / c(X^ , 0)dt; </°(t + r, X r °*)dr 



(2.12) 



for k — 1, . . . , n — 1 under [I]. The proof is almost the same as Theorem 5.1.9 in [T3] . 
Now we are prepared to state our main result whose proof is given in Section 15.31 

Theorem 3. Assume [A] -[I]. There are positive constants C n and rh n which are independent 
of e such that 



n-l 



u £ (t,x) - (u°(t,x) + JVt;°(t,x)) 



k=l 



< C n (l + \x\ lmn )e n , (t, x) E [0, T] x D. 



3 Semi-Group Representation 

In this section we construct a semi-group corresponding to (X^ ,x ) t and D, and give another 
form of (I2.12|) . We always assume [A]- [I] (or the generalized assumptions in Section I57T1 ) We 
only consider the case where c(t, x, 0) is non-negative and independent of t; we simply denote 
c(x, 0) = c(x). Let C®(D) be the set of bounded continuous functions / : D — > K such that 
f(x) = on dD. Obviously, Cj}(D) equipped with the sup-norm becomes a Banach space. 



For t E [0, T] and / E C$(D), we define P t u f : D 



P t D f(*) = E 



by 



■ /q c(X»>*)dv ) f(X^)l {TD(x0 ,^ t} 



We notice that Pff(x) is equal to u°(T — t,x) with the payoff function /. 



Theorem 4. The mapping P t : C^(D) 
contraction semi-group. 



Cjj(D) is well-defined and (P t )o<t<T is a 
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Proof. Let / G C° b {D). The relations P°f = f, P t D f\ dD = and sup \P t D f\ < sup |/| are 



D 



D 



obvious. The continuity of Pff is by Lemma 4.3 in [T7]. The semi-group property is verified 
by a straightforward calculation. ■ 



Note that (P t )t also becomes a semi-group on the set C p (D) of continuous functions /, 
each of which has a polynomial growth rate and satisfies f(x) = on dD. 
Let g E n m ' a > p . Observe that 



(T-t)AT D (X°<*) 



exp ( - / c(X°' x )dv J s(t + r, X r °' x )dr 



T-t 



exp ( - / c{X° v ' x )dv ) g(t + r, X^ x )l {TD{x o, x) > r} dr, 



0,x\ 



and we obtain 



E 



(T-t)Ar D (X°' x ) 









T-t 
T-t 



exp - / c W*)^ + r, X r °' x )dr 



E 



exp - / c(X^)dv g(t + r, X° r <*)l {TD{x o, x) > 



r} 



dr 



P r v g(t + r,-)(x)dr. 



Thus, under the assumption [H], we see 



P r D ^u°(T -t + r,-)(x)dr 

1 P r D J? l °P t D J(x)dr= f P t D r ^ l °P 7 D f(x)dr. 
o Jo 



Similarly we get the following. 
Theorem 5. For each k = 1, . . . , n — 1 

v° k (T-t,x) 



(3.1; 



E E ff ■ • ' f 1 p *-^^ 2 4° • • • pt^psm*i -dti. 

'=1 08 i )i = iCN , ,Eii8 1 =* 



(3.2) 



Proof. By (13. ip . we have the assertion for = 1. If the assertion holds for 1, . . . , k — 1, then 



!(T -*,*) = f P£{^V + E ^?-,«?}(T -t + to, -)(x)dt 
Jo i=i 



P t D t0 ^P t D J{x)dt, 
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k-1 I 

+EE E 

1=1 m=l ^ i )? 1 CN m ,E j ^=l 



JO J 



P t D - tn J? k °P t °f(x)dt 



-t ~^k to 

k 1 pt pti pti fh-\ 

'o Jo JO 



+EE E 

Pt D t 1 ^ 1 Ph-t 2 ^ 2 ^-t 3 -^'^ ' ' ' P^-tt'&piPt? f( x )dti -"dti 

= E E ff ' ■ ■ f " p ^4 ^-^ • • • tt^pfmdti 

Thus, our assertion is also true for fc. Then we complete the proof of Theorem [5] by mathe- 
matical induction. ■ 



In particular, when <i = 1, D = (/, oo), 0) = p, <r(x, 0) = cr and c(x) = with constants 
; G R and a > 
is well-known that 



/, /i G R and a > 0, the process Jf t °' x is explicitly represented as X^' x = x + pt + crS t , and it 



P(td(X°>*) > t|X°'*) = 1 - exp f- 2 ^ 0( f°' X S ) on {X t ^ > /}. 



Therefore, for g G Cp(D) we have 

/oo 
g(y)p(t,x,y)dy, (3.3) 

where 

. 1 2(g-0(y-Q (a-^-Mt) 2 

p(t, x, y) = (1 - e )e (3.4) 

V 2TT(J z t 

We remark that (I3.3P is useful for explicit evaluation of (13. ip . which is demonstrated in the 
next section. 



4 Application to Barrier Option Pricing under Stochas- 
tic Volatility 

Consider the following stochastic volatility model: 

dS £ t = (c - g)Sfdt + <7 t e S£<Z£ t \ S £ = S (4.1) 

da\ = e\(6 - a £ t )dt + svo\{pdB\ + y/\-pHB\), a e = a, 

where c,q>0,eE [0,1), A, 9, v > 0, p G [—1,1] and £> = (B X ,B 2 ) is a two dimensional 
Brownian motion. Here c and g represent a domestic interest rate and a foreign interest rate, 
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respectively when we consider the currency options. Clearly, applying Ito's formula, we have 
its logarithmic process: 

dX £ = (c-q- -{al) 2 )dt + a £ t dB], X £ = x = log S, (4.2) 



da\ = eX(9- a £ t )dt + eva\{pdB\ + y/l- p 2 dB 2 ), a £ = a. 
Also, its generator is expressed as 

= ( c - « - r 2 ) s + r 2 ^ + £ ^ 2 aS; + ^ " + 4"£- < 4 - 3 > 

In this case, J^f defined by (12. 5p is given as 

We will apply Theorem [12] to ( 14. ip with d = 2 and d' = 1 and give an approximation formula 
for a barrier option of which value is given under a risk-neutral probability measure as 

CS2L(T - t, e x ) = E [e^^f(S&)l {nL ^ ia .* )>T ^ } 

where / stands for a payoff function and L(< S) is a barrier price. 
u £ (t, x) = C|Ser( T - eX ) satisfies the following PDE: 

— + ^ £ -c)u £ (t,x) = 0, (t, x) G (0, T] x D, 
u £ (T,x) = f(x), xeD, { > 

u £ lt,l) = o, te[o,T]. 

where f(x) = max{e x — K, 0}, D = (/, oo) and I = logL. We obtain the 0-th order u° as 

u°(t, x) = P£J(x) = E[e- c ( T -*)/(X^)l {TD(x o„ )>r _, } ]. (4.6) 

Set a = c — q. Then P^_ t f{x) = C B ^ rrici (T — t, e x , a, a, L) is the price of the down-and-out 
barrier call option under the Black-Scholes model: 

l— 2a 2 

CiLUT - t, e x , a, a, L) = C BS (T - t, e x , a, a) - (j-j * C BS (t - t, L a, a) . (4.7) 

Here, we recall that the price of the plain vanilla option under the Black-Scholes model is 
given as 

C BS (T - t, e x , a, a) = e- q ^e x N(d 1 (T - t, x, a)) - e~ c(T ~ t] <KN(d 2 {T - t, x, a)), (4.8) 
where 

x - log K + at 1 n 

di{t,x,a) = -= h -aVt, 

ayjt 2 
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d 2 (t,x, a) = di(t, x, a) — a\fi 
N(x) = / n(y)dy, 



v In 



tr 



We show the following main result in this section. 

Theorem 6. We obtain an approximation formula for the down-and-out barrier call option 
under the stochastic volatility model 

C££«(r, e") = Cf£ rier (T, e», a, a, L) + «;?((), *) + 0(e 2 ), (4.9) 



u°(0,x) = e~ cT I I (1 < -~ )< t?(s, y)dyds, (4.10) 



1 2(rc-i)(j/-0 (j/-^-(c»-^ CT 2 ) 3 ) 2 




'o v2na 2 s 



= e a{T - t) puae x n{d 1 {T - t, x, a))(-d 2 (T - t, x, a)) 
+ 2e a( - T -Vpva ( %- J " Ln{ Cl {T -t,x,a))VT -t 



L 



2a 



- ea{T ~ t] P ua \Y ) Ln(d(T -t,x,a))ci(T -t,x,a) 

x { (r - t, ^, a, a) { 1 + (, - log L) (l - ^ ) ^ - (* 

+ A(^-a)e a(r - i) e :c n(c/ 1 (T-t,x,a))v / T :r t 



logLje-^^iV^^-^x^))} 



2a 



A(0 - (t) ( — ) e a{T - t] Ln{ Cl {T - t, x, a))VT~^t 



4a /, e x \ I e x \ / m Z/ 



^-*>A(0 - a)^ f log ^ J f ^ ) ( T - f , ^, a, a ) , (4.11) 



and 



. . 21 — x — log K + at 1 r- 

ci(i,x, a) = -= V -crvt. 

o\Jt I 

Proof. By Theorem [T2] and the equality (13.1 j) , we see the expansion 

T-t 



C^ iei (T - t, e x ) = C|f rrier (T - t, e x , a,a,L) + s f P? ^P°_ t _J{x)ds + 0(e 2 ). (4.12) 

J o 
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The first-order approximation term v®(t,x) = / P s D «5f 1 °P^„ t „ s /(x)(is is given by 

Jo 

= f * E[e- cs J??P°_ t _J(X T + s )]ds 
Jo 

= [ * E[e^^e^-^P^J(X T ^)]d8 
Jo 

= P s D J??P°_ t _J(x)ds, 
Jo 

where P f D is defined by (13.31) with the density (13.41) . that is, 

P t D f(x) = J - e f(y)dy. 

Define $(t,x) as 

■&{t,x) = ^P°_ t f{e x ) 

= e< T ^pva 2 —CiS TTicr (T - t, e x , a, a, L) + e< T ^\(9 - a)— C|^ rier (T - t, e*, a, a, L). 

A straightforward calculation shows that the above fucntion agrees with the right-hand side 
of (14. lip . Then we get the assertion. ■ 

Remark that through numerical integrations with respect to time s and space y in (I4.10p . 
we easily obtain the first order approximation of the down-and-out option prices. 

Next, as a special case of ( 14. ip we consider the following stochastic volatility model with 
no drifts: 

dSI = a £ t S £ t dBl S £ = S>0 1 (4.13) 
da\ = evo\(pdB\ + p 2 dB 2 ), a £ = a > 0. 

where e e [0, 1), p e [-1, 1] and B = (B 1 , B 2 ) is a two dimensional Brownian motion. In this 
case, we can give a slightly simple approximation formula compared with Theorem El 
By Ito's formula, the following logarithmic model is obtained. 

dX £ t =- l -{a £ t fdt + a £ t dBl X £ = x = \ogS ) 
daf = eisa £ t (pdBl + y/l - p 2 dB 2 ), a £ = a, 

This model is regarded as a SABR model with /3 = 1 and known as the log-normal SABR (see 
[T2]). Again, the barrier option price is given by 



^B^'rier(^' O ~ E [/(^T ) l{min < u < T 5g>L}] , 

where / stands for a payoff function and L(< S) is a barrier price. 
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The differentiation operators Jzf £ , Jzf® and the PDE are same as (I4.3p - (j4.5p with c = q = 
and A = 0. Also, the barrier option price in the Black-Scholes model coincides with (14. 7p with 
no drift, that is, 

C™ lat {T,S) = C BS (T,S)-(?) C BS (t. L 



L) \ ' S 

where C BS (T, S) is the driftless Black-Scholes formula of the European call option given by 
C BS (T, S) = SN(dt(T,logS)) -KN(d 2 (T,\ogS)) 

with 

, . . , . . x — log K + oHll 

di{t,x) = di(t,x,0) = -=. , 

a\Jt 

d 2 (t,x) = d 2 (t, x, 0) = di(t, x) — ay/t. 

Then, we reach the following expansion formula which only needs 1-dimensional numerical 
integration. 

Theorem 7. c££-r(r, e *) = CELer(T, e x ) + ev°(0, x) + 0(e 2 ), where 

v°{0,x) = -]^Tvpa{e x n{d l {T 1 x))d 2 {T,x)+Ln{c 1 {T,x))c 1 {T,x)} 

upL(x-l)log(L/K) f T (T-syl 2 ( c 2 (T-s,L/K) + c 2 (s,L/e*)\, 

exp as, 



2ira J s 3 / 2 

log(L7e-g) + aH/2 flogy + aHftV 

Cl[t ' x) = 77t ' C2< ' v) = V *St ) ■ 

Proof. By Theorem [T2] and the equality (13. ip , we see that the expansion 

C£S«(r, e x ) = C^JT, e x ) + sv^O, x) + 0(£ 2 ) 

holds with 

vl&x) = [ T ' ' P°_ t _ r J??P r D f(x)dr. (4.15) 
Jo 

Then, we have the following proposition for an expression of v®(0, x). The proof is given in 
Section 15.51 

Proposition 1. 

«i°(M = f ^fl^^>) ~ ^E[(T-r D (X°-))^ 2 ^P^ TD(x o,, ) /(0l^(^ ) <T } ]. 
We remark that the expectation in the above equality can be represented as 

\ E[(T - r D (X^))^a 2 ^-P^ rj5(x o,, ) /(0l { . D (xo.)<r } ] 
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= / ^^vpv 2 ^P?- s fm(s,x-l)ds, (4.16) 
where h(s, x — I) is the density function of the first hitting time to I defined by 

ft(s ' :c -' ) = VlS5F exp ( S5 I' (4 ' 17) 

Now we evaluate 



Note that 



and 



(i, e x ) = e :c n(d 1 (t, (4.18) 



Then we have 

, <9 2 



z/pcr 



dxda ' [ cr-s/i 



= -upae x n(d 1 (t,x))d 2 (t,x) (4.20) 

and 

vpa 2 ^— { (^] C BS (t, ^] ) = upaLn( Cl (t,x)) Cl (t,x). (4.21) 



dxda I \ L 
Combining (jOBjl . (jOQ| and f l4T2Tj) . we get 

vpa 2 -^-P t D f(x) = iypa{e x n(d 1 {t,x))(-d 2 {t,x)) - Ln( Cl (t,x)) Cl (t,x)} . (4.22) 
Substituting lj£22j) into f l4TTSj) . we have 

u P a2 inr P t D f( l ) = vp°Ln(d x (t,l)){-d 2 {t,l)) - paLn( Cl (t,l)) Cl (t,l) 
dxda 

= vpaLn{d x {t, l))(-(di(t, I) + d 2 (t, /))) 



1 / {I -log K + la 2 t} 2 \ f -2(1- log K) 

vpa exp 1 



Thus we obtain 



~\ m - ^(X°-))^a 2 ^-P^ Tc(x o. ) /(/)l {rD( xo.)< n ] 
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T (T-s) 1 _ {'-'°«*+^-*)} f -2(1- log K) 



o 



2vr V vVT 



-(l-x) {c-)+(- 2 /^} 2 

-e ds 



r 2 e 3 

vpL{x - I) \og{L/K) f T {T-s) 1 / 2 ( c 2 (T-s,L/K) + c 2 (s,L/e 3 



2na J ,3/ 2 ™V{ ! 2 ' ) ^ 

(4.23) 

By Proposition [H f l4.16j) . f)4.22p and (14.23)) . we reach the assertion. ■ 

Finally, we show a simple numerical example of European down-and-out barrier call prices 
as an illustrative purpose. Denote u° = C^ rrici (T, S) and v® = t>i(0, logs'). Then we see 

c£S-(r, S) ~ u° + £ < 

We list the numerical examples below, where the numbers in the parentheses show the error 
rates (%) relative to the benchmark prices of C^£ icr (T, S); they are computed by Monte 
Carlo simulations with 100,000 time steps and 1,000,000 trials. We check the accuracy of 
our approximations by changing the model parameters. Case 1-6 show the results for the 
stochastic volatility model with drifts (14.11) . and case 7 shows the result for the lognormal 
SABR model (Q51) . 

Apparently, our approximation formula u° + ev® improves the accuracy for C Bar '^ ier (T, S), 

and it is observed that ev® accurately compensates for the difference between C Bar '^ ier (T, S) 
and C|^ rier (T, S), which confirms the validity of our method. 



1. 



S = 100, a = 0.15, c = 0.01, q = 0.0, eu = 0.2, p = -0.5, 
e\ = 0.00, 9 = 0.00, L = 95, T = 0.5, if = 100, 102, 105. 



Table 1: Down-and-Out Barrier Option 



Strike 


Benchmark 


Our Approximation (w u + evy) 


Barrier Black-Scholes (w u ) 


100 


3.468 


3.466 (-0.05%) 


3.495 (0.80%) 


102 


2.822 


2.822 (0.00%) 


2.866 (1.57%) 


105 


1.986 


1.986 (0.01%) 


2.052 (3.36%) 



2. 



S = 100, a = 0.15, c = 0.01, q = 0.0, ev = 0.35, p = -0.7, 
e\ = 0.00, 9 = 0.00, L = 95, T = 0.5, if = 100, 102, 105. 
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Table 2: Down-and-Out Barrier Option 



Strike 


Benchmark 


Our Approximation (w u + ev^) 


Barrier Black-Scholes (w u ) 


100 


3.421 


3.423 (0.07%) 


3.495 (2.18%) 


102 


2.753 


2.757 (0.18%) 


2.866 (4.13%) 


105 


1.885 


1.890 (0.23%) 


2.052 (8.88%) 



3. 

S = 100, a = 0.15, c = 0.05, q = 0.0, ev = 0.35, p = -0.7, 
eX = 0.00, 9 = 0.00, L = 95, T = 0.5, K = 100, 102, 105. 



Table 3: Down-and-Out Barrier Option 



Strike 


Benchmark 


Our Approximation (u u + ev^) 


Barrier Black-Scholes (n u ) 


100 


4.352 


4.349 (-0.07%) 


4.399 (1.06%) 


102 


3.585 


3.586 (0.02%) 


3.665 (2.24%) 


105 


2.560 


2.563 (0.11%) 


2.696 (5.31%) 



1. 

S = 100, a = 0.15, c = 0.05, q = 0.1,eu = 0.2, p = -0.5, 
e\ = 0.00, 9 = 0.00, L = 95, T = 0.5, K = 100, 102, 105. 



Table 4: Down-and-Out Barrier Option 



Strike 


Benchmark 


Our Approximation (w u + ev^) 


Barrier Black-Scholes (w u ) 


100 


2.231 


2.224 (-0.31%) 


2.268 (1.64%) 


102 


1.758 


1.754 (-0.27%) 


1.812 (3.02%) 


105 


1.172 


1.168 (-0.31%) 


1.243 (6.05%) 



5. 

S = 100, a = 0.15, c = 0.01, q = 0.0, ev = 0.2, p = -0.5, 
e\ = 0.2, 9 = 0.25, L = 95, T = 0.5, if = 100, 102, 105. 
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Table 5: Down-and-Out Barrier Option 



Strike 


Benchmark 


Our Approximation (w u + ev^) 


Barrier Black-Scholes (w u ) 


100 


3.523 


3.517 (-0.16%) 


3.495 (-0.77%) 


102 


2.891 


2.888 (-0.09%) 


2.866 (-0.85%) 


105 


2.066 


2.065 (-0.06%) 


2.052 (-0.64%) 



S = 100, a = 0.15, c = 0.01, q = 0.0, ev = 0.2, p = -0.5, 
£ \ = 0.5, 9 = 0.25, L = 95, T = 0.5, K = 100, 102, 105. 



Table 6: Down-and-Out Barrier Option 



Strike 


Benchmark 


Our Approximation (w u + ev^) 


Barrier Black-Scholes (w u ) 


100 


3.587 


3.594 (0.20%) 


3.495 (-2.55%) 


102 


2.976 


2.987 (0.39%) 


2.866 (-3.68%) 


105 


2.170 


2.183 (0.59%) 


2.052 (-5.41%) 



S = 100, a = 0.15, c 
£ \ = 0.0, 6 = 0.0, L -- 



-- 0.0, q = 0.0, ev = 0.2, p = -0.5, 
95, T = 0.5, K = 100, 102, 105. 



Table 7: Down-and-Out Barrier Option 



Strike 


Benchmark 


Our Approximation (w u + ev^) 


Barrier Black-Scholes (w u ) 


100 


3.261 


3.258 (-0.09%) 


3.290 (0.90%) 


102 


2.640 


2.639 (-0.02%) 


2.686 (1.78%) 


105 


1.841 


1.841 (0.01%) 


1.911 (3.77%) 



5 Appendix 

5.1 Generalization of Main Results 

There are several cases in practice that our assumptions [A]-[B] are not satisfied. Hence, 
in this section we weaken the assumptions. Let d! G {!,..., d}, and we regard X\ ,x ' 1 as 
logarithm of the underlying asset prices for i < d', and as parameter processes such as those 
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for a stochastic volatility and a stochastic interest rate for i > d! . For a technical reason 
introduced later, we assume I C [0, oo) in this section. 

Let (Q, J 7 , (J-j)t, P) be a filtered space equipped with a standard Brownian motion (B t )t- 

Set 

gi (l/je) = \y l \b\<y),e) + \j^y\v{y),e))^ i<d', 
{ 6*(7r(2/),e), i > d', 

ly ' j " I ^(7r(i/) >e ), ixf, 
where n(y) = (log y 1 , . . . , log /, y d '+\ . . . , y d ) G R d . 

[A'] For each e G 7 it holds that cr 4J (-, e), b l (-, e) G £ and that a tj (-,e), b\-,e) and c(7r(-),e) 
are also in C, that is, they are extended to be locally Lipschitz functions (with respect 
to the parabolic distance) defined on M. d . Moreover, there exists a solution (X £ ' x ) t of 
SDE (ll.3p and for any m > there are m', C > such that 

sup E[|y^H < cr- x (i + m w ), g [o,t] x [o,oo) d ' x R d ~ d ', e e /, (5.1) 

0<r<t 

where Y t £ ' y = t(X £ ' n(y) ) and l(x) = (e x \ . . . , e xd \ x d ' + \ . . . , x d ) G M d . 

[B'] The function /(x) is represented by the continuous function / : M. d — > R as /(x) = 
There exists Cf > such that \f(y)\ 2 < C;(l + |y| 2m ), y G M d . Moreover, 
f(x) = 0onR d \D. 

[C] In addition to the condition [C], there is a constant A £ A > such that c(t, < 
^(l + |x| 2m ). 

Note that Ito's formula implies that (Y t £,y ) t is a solution of 

dY t £ ' v = b(Y £ ' y , e)dt + a(Y £ ' y , s)dB t , 

Y £ ' y - 7/ 

Although Theorem 3.1 in [T7] no longer works under [A']-[B'], we can charaterize u £ (t,x) as 
the solution of (12.31) in the viscosity sense. To see this, set 

i,j — 1 i — 1 

d 

where d* J = a lk a^ k . Moreover, define 

k=l 

D = {y G R d ; y { > 0, i = 1, . . . , d' and ir{y) G D} 
and u e (t,y) = n £ (t,vr(y)) ((t,y) G [0,T] x D), ((t,y) G [0,T] x dD). 
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Theorem 8. Assume [A']-[B']. Then, u e (t,x) is a viscosity solution of (I2.3p . Moreover, 
u £ (t,y) is a viscosity solution of 

-^u £ {t, y) - J? £ u £ (t, y) = 0, (t, y) G [0, T) x D, 

H e (T,y)=f(y), xeD, (5 ' 2) 



satisfying 



u £ (t,y) = 0, (t,y)e[0,T]xdD 



sup _\u s (t,y)\/(l + \y\ 2m ')<oo. (5.3) 

(i,y)g[0,T]x£> 



Proof. The latter assertion is by the similar argument to the proof of Proposition Then, 
the simple calculation gives the former assertion. ■ 

Applying Theorem 8.2 in [I] and Theorem 7.7.2 in [15] for (15.21) . we have the following 
theorem. 

Theorem 9. Assume [A']-[C'] and [D]-[E]. Ifw e (t,y) is a viscosity solution of (15.21) satis- 
fying the growth condition (15.31) . then u e = w e . 

Let U m > a <P be the same as U m ' a ' p replacing O with 

\g(t,x)\ < M 9 (t) {1 + \t(x)\ 2m } , t G [0,T), x,y G D. 

Moreover we define Q m ' a ' p similarly to Q m > a ' p , replacing H mA2 and U m > a ' p in the definition 
with n m ' ' 2 and n m ' a ' p , respectively. 

[H'] The condition [H] holds replacing Q m ^ p with Q m ^ p '. 

[I'] The condition [I] holds replacing G m ' a ' p with G m > a ' p . 

The following theorem gives a generalization of Theorem [3j The proof is in Section 15.41 

Theorem 10. Assume [A']-[C], [D]-[G] and [H']-[I']. Then there are positive constants C n 
and m n which are independent of e such that 



n-l 



u e (t,x) - (u°(t,x) + J2^l(t,x)) 



< C n (l + \L{x)\ 2 ^)e n , (t,x) G [0,T] x D. 



k=i 

Here we give another version of generalized assumptions. 

[D'] The domain D is given as D = Ux R d ~ d ', where U is a domain in M. d whose boundary 
dU satisfies the outside strong shpere property. 

[E'] The condition [E] holds for e ^ 0. Moreover a ij (x, 0) = V[x, 0) = for i = d! + 1, . . . , d, 
j = 1, . . . , d and for each compact set K C D there is a positive constant fio,K such that 

d! 

i^, K \t,\ 2 < a '^ x ' °) - ^l^l 2 for x e K and £ e Rd - 
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[F'] For each y E R d d ' , the inequality 

maxsup{|<7^((x,y),0)| + |6*((x, y), 0)| + |c((x, y), 0)|} < oo 

holds, where (x, y) = (x 1 , . . . , x d \ y 1 , . . . , y d ~ d ') E R d and there exist A 3 (y) > and 
a E (0, 1] such that 

|<7*'((ar, y), 0) - ^'((x', y), 0)| + |6 4 ((a;, y), 0) - b\(x', y), 0)| + |c((x, y), 0) - c((a/, y), ( 
< A 3 (y)|x-x'r, G [0,T] x U,y E R d ~ d ' . 

[H"] The condition [H] holds replacing g m > a ' p with Q ' ' , where 



<T W ' = {yGC 1 ' 2 ([0,T)x J D)nC([0,T] x 5) 



and is the same as U m ' a <P replacing D C M d in the definition with £/ C R d ' . 

Here h(-,-;y) denotes the function [0,T] x U 3 (t,x) i — >■ h(t,(x,y)) G R for /i = 
dg /dx l , d 2 g /dx l dx j . 

[I"] The condition [I] holds replacing g m ' a ' p with ZT' Q,P . 

Theorem [2] implies the next theorem. 

Theorem 11. Assume [A']-[F'] and [G]. Let g E T-L^' a,p for some p > 1/a. Then for each 
fixed y = {x % ) d =d , +1 , the following PDE 

^v(t, x) + Stfvit, x) + g(t, x- y) = 0, (t, x) E [0, T) x U, 

v(T,x) = 0, xEU, 
v(t,x) = 0, (t, x) G [0, T] x dU 

has a classical solution v satisfying (12. lip , where 

i,j=l i=l 

Moreover, if w is another classical solution of f)2.10p which satisfies \w(t,x)\ < C exp((3\x\ 2 ), 
(t, x) E [0, T] x D, for some C, (3 > 0, then v = w. 

Using the above theorem instead of Theorem [2] itself, we can prove the following theorem 
similarly to Theorem [TU1 

Theorem 12. Assume [A']-[F'], [G] and [H"]-[I"]. Then, the same assertion of Theorem 
\M holds. 
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5.2 Proof of Theorem H 

We consider the following PDE which is equivalent to (12. 10p with changing variable t to 
T-t 

d 

-— v(t, x) + Sf°v(t 7 x) + g(t, x) = 0, (t, x) G (0, T] x D, 

v(Q,x) = 0, xeD, ( 5 - 4 ) 

v(t,x)=0, (t, x) G [0, T] x 3D. 

We define H m ' Q ' p as the same as U m ^ v replacing [0,T) in the definition with (0,T]. 
We divide the proof of Theorem [2] into the following two propositions. 

Proposition 2. For any g, a classical solution of (15.41) is unique in the following sense: if v 
andw are classical solutions of (15.41) and \v(t,x)\ + \w(t,x)\ < C exp((3\x\ 2 ) for some C, [3 > 0, 
then v = w. 

Proposition [2] is obtained by the same argument as the proof of Theorem 2.4.9 in [5]. 

Proposition 3. There exists a classical solution v of (15.41) for g G "H m,a ' p with p > 1/a. 
Moreover, (12. lip holds. 

Proof. By Levi's parametrix method, we can construct the fundamental solution T(t,x;T, £) 
for the operator L = —d/dt + Jzf°, that is, 

ft 



W g (t,x)= f [ r(t,x;T,£)g(T,Z)dZdT 

JO J D 



is continuous in (t, x), continuously different iable in x for g G C([0, T] x D). When g is Holder 
continuous in x uniformly in t G [0, T], then we see that W g is a solution of (15.41) (See Theorem 
1.5.8-1.6.10 in [5]. For more details, please refer to Chapter 1, Section 2-6.) However functions 
in % m ' a ' p may not have the regularity at t = 0. So we generalize the argument in Chapter 1 
of [5]. We remark that T(t, x; r, £) is given by 

r(t,x;<r,£) = Z(t,x]T,£) + / / Z(t,x;a,r))$((T,r);T,£)dr)d(T, 



D 



where 



Z( t , x; r, £) = exp ( ^ " W " ? 



(4vr(t - r)) d / 2 4(t- r) 

and s; r, £) is the solution of 

$(t,x;r,£) = LZ(t,x;T,£) + J J LZ(t,x;a,r])®(o',r];T,£)dr]do'. 

Fix any # G H m ' a ' p . We can divide as W g = V g + U g , where 

Vg(t,x)= / / Z(t,x;T,£,)g(r,£)d£dT, 

JO J D 
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U g {t,x)= [ [ Z(t,x]T,g)g(T,£)d£dT, g(t,x)= [ [ $(t,x]T,£)g(T,£)d£dT. 

JO J D JO Jd 

We remark that V g , U g and g are well-defined by virtue of (4.9) and (4.15) in j5] and the 
property of g. Take /3 G (a — 1/p, a). By Theorem 1.4.8 in [5], we see that 

* f , x- r, - $(t, y; r, 1 < t; J^lL^v, exp ?L + exp ' 



( t _ r )(d+2-(a-/3))/2 |^ ^ t _ r y ' t _ r 

for some C, A > 0. Hence, 

\g(t, x) - g(t, y) \ < C f —J^L—dr\x - yf 



'O (* - T)M«-/9)/2 
Vp / /-t \ 1/9 

< a ' 



y (AF(r)) p dr) (j(t-T)-^ a -^ q drj \x-yf, t G (0, T], x,y E D 

for some C" > by virtue of the Holder inequality, where q > 1 is given by 1/p+l/q = 1. Since 
(1 — (a — 0)/2)q is smaller than 1, we see that g(t,x) is /3-H61der continuous in x uniformly 
in £ G (0,T]. Then, Theorem 1.3.3-1.3.6 and the equality (4.2) in Chapter 1 of [5] imply that 
Ug(t,x) G C 1,2 ((0,T] x D) and 



LU g (t,x) = -g(t,x) 



f I UX>^°)-^ )) 

Jo Jd [ z ij=1 



2 



dx % dxi 



-b\x, 0) A - c(x, 0) }>Z(t, x; r, £)£(r, £)d£dr 



= -£(*,z)4- / / LZ(t,x]T,£)g(T,£)d£dT. (5.5) 
Jo JO 

For the volume potential V g , we follow the proof of Theorem 1.3.4 in [5] to find that for 
any compact set K C D 



E 



a 2 



dx l dxi 
for some > 0, where 



^~77 te (°' T ]' re(0,t), fie (1- a/2,i; 



J g (t,x,r)= / Z(t,x;T,£)g(r,£)d£. 
Jd 



Hence, the dominated convergence theorem implies that V^(t, x) is twice continuously differ- 
ent iable in x. Similarly, we get V g G C 1,2 ((0,T] x D) and 

LV g (t,x) = -g(t,x)+ [ [ LZ{t,x;T,£)g{T,£)d£dT. (5.6) 

Jo Jd 

Combining fl5.5p - fl5.6p . we obtain 

LW g {t, x) = -g(t, x) - [ [ {$(*, x- t, - LZ(t, x; r, £)} g(r, £)d£dr 

Jo Jd 
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ft 

+ 



/ / / / LZ(t,x;T,^)g(a,rj)drjdad^dT 
Jo J d Jo J D 



-g(t,x)- / {$(t,x;T,Z)-LZ(t,x;T,£) 
Jo Jd 



J J LZ(t,x;a,r])^(a,T];T,()dT]da\g(T,^)d(dT 



-g{t,x) 



which implies that W g is a solution of (15.41) . Moreover, since g G 'H m ' a ' p , using the inequality 
(6.12) in p.24 of [5j, we get 

< C" f f T (M 9 (r)) p dT) (1 + |x| 2m ) 



for some C",C", \' > 0. Then, we complete the proof of Proposition [3J 

5.3 Proof of Theorem G2 

First, we generalize the definitions of Jzf fc °, and We define 

i f i _ d _ rl f)k n ij f)2 

v i s j — 1 



n— 1 n—2 ( n— 1 



(5.7) 



fe=l fe=l k Z=fc+1 

+e"- 1 jSf>£_ 1 (*,*). 
We consider the following Cauchy-Dirichlet problem: 

- |U(t, x) - J? £ v(t, x) - g £ n (t, x) = 0, (t, x) G [0, T) x D, 

t;(T,x)=0, x £ D, 

v(t,x)=0, (t, x) G [0, T] x 

For £ 7^ 0, we define v e n = [u 6 — {u° + ^ £ fc f °(t, x)}]/5™. Obviously, we see 

k=l 
n-l 

u £ {t,x) = u\t,x) + ^e k vl(t,x) + e n v £ n (t,x). (5.8) 

k=l 
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Proposition 4. The function v e n is a solution of (15. 7p . 

Proof. It is obvious that v%[T,x) = for x G -D and v £ n (t,x) = for (t,x) G [0,T] x 
Apply Taylor's theorem to (12. 3p to ovserve 



n-l 



J^ e u £ (t, x) = < ^° + £k ^k + \ u£ (t> x). 



k=l 



Since u° is the solution of ( 12. 3 p with e = 0, we get 

5 



^-u°(t,x)+Sf°u°(t,x) = 0. 



Similarly, by Theorem El we have 




k-1 



dt 



(t, x) + Sf°v° k (t, x) + JSf°u°(*, x) + J2 -Zk-Piit, x) = 0. 



i=i 



Combining (j5.8j) — (j5.11j) and Theorem [TJ we obtain 







n-l 



e n { -<(t, x) + J2f°<(t, x) + ^>°(t, x) + £ *) 



2n-2 



n-l 



+ Y, £k \^-n<(t,x)+j?y k _ n (t,x)+ Yl *S-iAt, 



x 



k=n+l 
2n-l I rZO „.e, 



_2n c2i£ n ,£ , 



and thus, 



+e*- 1 JCXft x) + J?*vU(t, x) + e*J2^<(t, a) = 



^<(t, x) + .i? £ <(t, x) + <£(*, x) = 0. 



c(X^,e)dv)g E n (r + t,X^)dr 



(5.9) 



(5.10) 



(5.11) 



This implies the assertion. 
Set 

" nr D {X e ^)/\(T-t) 

v n {t,x)=E / exp 



By [F]-[I], we find that there are C n > 0, m n G N which are independent of e and the function 
M n G C([0, T)) n L p ([0, T), dt) determined by it , u°, . . . , v$_ t such that 



|<£(t,x)| <C n M n (t)(l + |x| 2 ^). 
The inequalities (12. 2p and (I5.12p imply 



K(t,x)\ < C' n J M n {r)dr{\ + |x| 2 ™") 
for some C > which is also independent of e. 



(5.12) 



(5.13) 
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/ = v e . 

n n 



Proposition 5. v\ 

Proof. The assertion is easily obtained by the similar argument to Theorem 5.1.9 in 



n-l 



Proof of Theorem^ By (15. 8 p and Proposition |5j we have u £ (t, x) — (u°(t, x) + e k v^{t, x)) 



k=l 



e n v^(t,x). Our assertion is now immediately obtained by the inequality (I5.13p . 



5.4 Proof of Theorem (TH 



Let v e n and v e n be as in Section 15.31 Thanks for the assumption I C [0,oo), the same 
argument as the proof of Proposition|4]tells us that v e n is a viscosity solution of (15. 7p . Moreover, 
we have the next proposition. 

Proposition 6. The function v e n is a viscosity solution of (15. 7p . 

Proof. Until the end of the proof we suppress e in the notation. First, we check the continuity. 
By the similar argument to the proof of Lemma 4.2 in [17], we see that v n is continuous on 
[0, T) x D. Moreover, by ( 15TT3|) . we get 



\2m\ 



sup \v n (t, x) | < C n (l + sup \x 
xeKnD xeK {Jo 



T ,t 

M n (r)dr- / M n (r)dr 



0, t^T 



for any compact set K C M. d . Thus, v n is continuous on [0, T] x D. 

Next, we show that u n is a viscosity subsolution of (15. 7p . Take any (t,x) G [0,T) x D and 
let </? be C 1,2 -function such that v n — ip has a maximum at (t,x). We may assume that <p 
and its derivatives have polynomial growth rates in x uniformly in t. By the Markov property, 
we have 



E[J(h)v n (t + h,X x h )l 



{T D {X*)>h} 



]=E 



(T-t)A(T D (X* +h )+h) 



J{r)g n (t + r,X*)drl {TD ( X *)>h} 



where J(r) = exp f- jf c(X*,£)d^. Since T D (X x +h ) = t d (X x ) - h on {r D (X x ') > /i}, we 



obtain 

E [J(/i)v„(t + h, Xl)l {TD{Xx) > h} ] = v n (t, x)-E 
where 



J(r)g n (t + r,X x )drl 



A h = E 



(T-t)AT D (X*) 



J(r)g n (t + r 7 X x )drl 

JO 

Therefore, 

ip(t, x) = v n (t, x) = E [J(h)v n (t + h, Xl)l {TD{x -)>h}\ + E[ / g n (t + r, X*)drl {TD{x *)>h}) + A h 

./o 
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rh 

< E [J(h)<f(t + h, X%)l {TDiX ')>h}\ + E[ / Quit + r, X?)drl {TD(x *)>h}} + A h 

Jo 



Applying Ito's formula to J(r)(p(t + r,X*), we get 



I 

h ., 



E 



| ^ + ^ + r,X*) + g n (t + r,X*)} l {rD (x-)>/.} 



< A h -ip(t,x)P(T D (X*)<h) ^ (514) 

By f)5.12p and the Schwarz inequality, we have 

\A h \ < C;'(l + \i{x)\ 2m ) f M T dtP{r D {X x ) < h) 1 ' 2 

Jo 

for some C" > 0. Using (15. ip and the Chebyshev inequality, we obtain 

/"(/// 

P(r D (X :c ) < /i) < E[ sup \X X -x\> dist(x, < -r— r J \ 7 ^ E[ sup |X r a - x| 8 ] 

0<r<h dlSt(X, (XL*) 5 r£[0,/i] 

(-11111 

dist(x,(9D) 8 

for some C"', C^", m > 0. Thus, letting h — > in (I5.14p . we see that 

d 

-— tp(t, x) - 3?<p(t, x) - g n (t, x) < 0. 

Hence, v n is a viscosity subsolution of ( 15. 7p . By the similar argument, we also find that v n is 
a viscosity supersolution. By the definition of v n , we easily get v n (T,x) = for x G D and 
v n (t, x)=0 for (t, x) E [0, T] x <9£>. ■ 

To see the equivalence v e n = v £ n , we need to give a new proof of Proposition \5\ under the 
assumptions of Theorem [TU1 

n-l 

Proof of Proposition^ Set u 6 n (t,x) = u°(t,x) + e k v^{t, x) + e n v £ l (t, x). The analogous 

fc=i 

argument of the proof of Proposition |4] implies that u e n is a viscosity solutions of (12. 3p . We 
easily see that u £ n has a polynomial growth rate in x uniformly in t. Then, Theorem [9] leads 
us to u e n = u £ . This equality and (I5.8P imply the assertion. ■ 

Now, we obtain the assertion of Theorem [10] by the same way as that of Theorem [3j 

5.5 Proof of Proposition [I] 

First, we notice the following relation: 

J??P t D f(x) = upoH - ^ P t D f(x). (5.15) 
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Then, using the relations ^ P t D f (x) = ^ P t D f (x) and 

d 



if P?_J{x) = 0, 



at 



we get 



Also, we have 



d_ 

Of 



T-t 







T — t ~ D 



2 



■^i P?-tf i x ) 



(5.16) 



+ if / P^_ t _ r ^pa 2 — -P r D / (x)dr = -J2?P£.J(x), x G (Z, oo). (5.17) 



Therefore, the function 

r-T-t 



T)(t,x) 







dxdcr 



d 2 



dxdu 



T-t ~ 



nfl 



(5.1* 



satisfies the following PDE 



/ ^ + J2f° ) 77(t,a;) = 0, 



77 (T, x) = 



7/(M) 



■?—±j??p°_ t f(l), te[0,T). 



(t,x) e [0,T) x (Z,oo), 
x E [I, oo), 



Then Theorem 6.5.2 in [6] implies 



1 B 2 

'~^0,x\\ .2 a pD 



^(0,x) = --E[(T-r D (X^))^ 



By (I5.18P and ( I5.19p . we get the assertion. 



}]• (5-19) 
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